Quantum sensitivity limit of a Sagnac hybrid interferometer based on slow-light 

propagation in ultra-cold gases 
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The light-matter-wave Sagnac interferometer based on ultra-slow light proposed recently in 
(Phys. Rev. Lett. 92, 253201 (2004)) is analyzed in detail. In particular the effect of confining 
potentials is examined and it is shown that the ultra-slow light attains a rotational phase shift 
equivalent to that of a matter wave, if and only if the coherence transfer from light to atoms as- 
sociated with slow light is associated with a momentum transfer and if an ultra-cold gas in a ring 
trap is used. The quantum sensitivity limit of the Sagnac interferometer is determined and the 
minimum detectable rotation rate calculated. It is shown that the slow-light interferometer allows 
for a significantly higher signal-to-noise ratio as possible in current matter-wave gyroscopes. 

PACS numbers: 42.50.Gy, 03.75.-b, 42.81.Pa 



I. INTRODUCTION 

In contrast to inertial motion, rotation of an object is 
absolute in the sense that it can be defined intrinsically, 
i.e. independent of any inertial frame of reference. Rota- 
tion can be detected e.g. by means of the Sagnac effect 
0, i.e. the relative phase shift A</> rot of counterpropa- 
gating waves in a ring interferometer of area A attached 
to the laboratory frame rotating with angular velocity ft. 
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where A is the wavelength and v the phase velocity of 
the wave. Depending on the nature of the wave phe- 
nomena employed, one distinguishes two basic types of 
Sagnac interferometer: laser 0, 0, 0] and matter-wave 
gyroscopes The Sagnac phase shift per unit area in 
a matter-wave device exceeds that of laser-based gyro- 
scopes by the ratio of rest energy per particle to pho- 
ton energy mc 2 jhw which for alkali atoms and optical 
photons is of the order of 10 11 [E0- Despite this very 
large number, matter-wave gyroscopes have only recently 
reached the short-time sensitivities of laser based devices 
This has mainly two reasons: First of all, laser- 
based gyroscopes, especially fiber-optics interferometer, 
can have a much larger area than matter-wave systems 
[Tof . Secondly the large flux of photons achievable in op- 
tical systems leads to a much lower shot-noise level as 
compared to matter- wave devices 0, . Thus in order 
to make full use of the much larger rotational sensitivity 
per unit area in a matter-wave device one needs to find 
ways to increase (i) the interferometer area and (ii) the 
particle flux. While a substantial increase of the inter- 
ferometer area in matter- wave devices is difficult, the use 
of novel cooling techniques has lead to high-flux atom 
sources which improved the performance of atom inter- 
ferometers 0, ^| . With particle throughputs which can 
now reach 10 8 s _1 as compared to a few atoms per second 
in the first atomic interferometers, the noise level is how- 



ever still much higher than that achievable in fiber op- 
tics gyroscopes with photon counting rates on the order 
of 10 16 s _1 Continuously loaded Bose- Einstein 

condensates (BEC) could provide a source for coherent 
atoms with larger flux values, and substantial progress 
has been made in this direction over the past few years 

m 

We recently proposed a light-matter-wave hybrid in- 
terferometer based on slow-light propagation in ultra- 
cold gases of 3- level atoms |14|. We argued that this 
interferometer would combine the large rotational phase 
shift of matter-wave systems with the large area typi- 
cal for optical gyroscopes. To this end the simultaneous 
coherence and momentum transfer associated with ultra- 
slow light in cold atomic gases with electromagnetically 
induced transparency (EIT) E3 was utilized. As the re- 
duction of the group velocity of light in 3- level EIT media 
is based on the change of character of the dressed eigen- 
modes of the systems from electromagnetic excitations 
to atomic Raman excitations |l6| , light waves can coher- 
ently be transformed into matter waves. These matter 
waves pick up a Sagnac phase shift per unit area which is 
orders of magnitude larger than the corresponding value 
for electromagnetic fields. 

In the present paper we present a detailed theoreti- 
cal description of the light-matter-wave hybrid interfer- 
ometer. In particular we discuss the effect of confining 
potentials for the atoms. We find that in contrast to 
the case of an infinitely extended medium or of periodic 
boundary conditions, which have been assumed in |l4j . 
the wavefunctions of all three internal states acquire the 
same matter- wave contribution to the Sagnac phase when 
in motional equilibrium with a trapping potential |17| . 
As a consequence the matter-wave contribution to the 
rotational phase shift vanishes. Only if periodic bound- 
ary conditions for the ground-state wavefunction can be 
maintained a nonvanishing matter-wave contribution to 
the rotational phase shift emerges. This can be realized 
e.g. in a circular-waveg uide BEC El El. The need for 
a circular atomic waveguide puts more stringent restric- 
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tions to the possible interferometer area then assumed 
in Q and thus partially invalidates the advantages of 
the hybrid interferometer stated in that paper. We will 
show however, that despite this restriction the minimum 
detectable rotation rate at the shot-noise limit can ex- 
ceed the current state of the art. It corresponds to that 
of a matter-wave gyroscope with a rather large particle 
flux given by the high density of the ultra-cold gas, e.g. 
a BEC, multiplied by the recoil velocity. To determine 
the quantum sensitivity limit of the hybrid interferometer 
the saturation of the Sagnac phase shift with the probe- 
light intensity as well as probe-field absorption will be 
taken into account. It will be shown that the Sagnac 
phase attains a maximum value for a certain value of 
the probe power. Optimum parameter values for a maxi- 
mum signal-to-noise ratio (SNR) will be determined and 
the minimum detectable rotation rate fi m i n per unit area 
derived. 



II. DYNAMICS IN THE ROTATING FRAME 

An intrinsic sensor attached to the laboratory frame 
detects the rotation of the frame without any reference 
to some other, non-rotating frame of reference. Thus it 
is most natural to describe this system from the point 
of view of a co-rotating observer. We will give here a 
microscopic description of the gyroscope which consists 
of an ensemble of three-level atoms with internal states 
|1), |2), and |3) in a ring interferometer, coupled by two 
laser fields with (complex) Rabi frequencies il c and f2 p 
in a Raman configuration as shown in Fig. ^ (left). The 
probe field fl p is assumed to propagate clockwise and 
counter-clockwise with respect to the rotation axis e z 
with its beam path bound to a circle of radius R as 
depicted in Fig. \I\(right). The control field f2 c , which 
is assumed to have a much larger Rabi frequency than 
the probe field, propagates in a different direction such 
that the corresponding wavevectors are (nearly) perpen- 
dicular. The ensemble as well as the laser sources are 
attached to the laboratory frame rotating with angular 
velocity O(t) = Q(t)e z [2(|. The center-of-mass motion 
of the atoms shall also be confined to the periphery of the 
circular loop. Furthermore, it is assumed that <C c 
such that non-relativistic quantum mechanics applies. 

Under conditions of two-photon resonance, the control 
field f2 c generates EIT for the probe field associated with 
a substantial reduction of its group velocity 0, |2^, |23j . 
The group velocity reduction which is due to the cou- 
pling of the probe light to the atomic Raman coherence 
corresponds in a quasi-particle picture to the formation 
of so-called dark-state polaritons, a superposition of light 
and matter degrees of freedom |24( • The smaller the 
group velocity the larger the contribution of the matter 
component in the polariton. 

The atoms are here described in second quantization 
by three Schrodinger fields 4'i(r, t), ^(r, t), and ^(r, t) 
corresponding to the three internal states. In order to de- 
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FIG. 1: (left) atomic level scheme, p denotes the momentum 
of the atoms and k p the wavevector of the probe-field along 
the periphery (x) of the circular loop depicted in the right 
part of the figure, fcjj in the component of the control- field 
wavevector along x, and Afc = k p — k\. (right) Schematic set- 
up of the hybrid Sagnac interferometer with vapor cell (grey 
box) attached to the rotating frame with angular velocity Q. 



scribe the propagation of the probe light and the three 
matter-wave fields in the co-rotating frame, we need to 
transform the Hamiltonian of the system into the rotat- 
ing frame. 

As the starting point we choose the standard atom- 
light interaction Hamiltonian of quantum optics in 
Coulomb g^uge and after the Power- Zienau-Wolley trans- 
formation [25l |. Adding the free Hamiltonian of a 3- 
component non-relativistic Schrodinger field, the Hamil- 
tonian reads in a non-rotating frame 

H = + h(f) + h(D 

h 2 



= ^|d 3 r*t(r) 



ii 



V^+^ + K ext (r,*) 



H(r) 

eo 



(V x A ± (r 



(2) 



+ J2 [ d 3 r %(r) [cV • (ri(r) - E cxt (r, t)))] tt tf (r). 

i?( A ) describes the motion of atoms in an external, 
possibly state- and time-dependent trapping potential 
V^ xt (r, t), Efj, — fthj^ is the energy of atoms in state 
The free Hamiltonian of the radiation field is denoted 
by , where Ai (r) is the transverse part of the vec- 
tor potential and II(r) its conjugate momentum. Finally 
describes the interaction of the atoms with the quan- 
tized electromagnetic field and additional external fields 
in dipole approximation, where is the vectorial dipole 
matrix element between internal states and \v). For 
notational simplicity we will drop the subscript "_L" in 
the following. 

The transition to a frame rotating with angular veloc- 
ity tt(t) is done via the unitary transformation 



U(t)=cxp[-- 



dr SI(t) ■ L 



(3) 



where L is the total angular momentum operator of light 
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and matter. By choosing f2 — Q e 2 we restrict ourselves 
to a rotation about the fixed z-axis. In this case only the 
vector component parallel to that axis, i. e. 



and 



oE d3r [fti(Mi) + (Mi)nj 



(4) 



is relevant. In eq. |@J the index /x denotes the three inter- 
nal states and the index j the three spatial dimensions. 
The Hamiltonian in the rotating frame is hence given by 



H mt = U{t)HU\t) + n(t)L z 



(5) 



Since L z an d £z commute, the unitary transforma- 
tion eq. © can be decomposed into two operators which 
act on the matter-wave and on the electromagnetic field 
respectively. One finds 



+E/ dV *K r ') 



2m 



/r v' 2 + hu„ + y; xt (r') 



(6) 
Mr'). 



(7) 



J2 / dV *t (r') [<!„„ • (n(r') - E cxt (r',i)) 1 *„(r') 



Here the prime denotes that the variables are given with 
respect to the rotating coordinates 



f dTe v Rn(T), 



(8) 



with R being the distance from the rotation axis. For all 
field operators F € {^n, A} holds: 



J7F(r)C/t = F(r + / dr e v i?^(r) 

V J to 



(9) 



The center-of-mass dynamics of the matter-wave fields 
is then governed by the following Heisenberg equations 
of motion in the co-rotating frame 



ih(d t + n(t)9 v )^(r',t) = 



(10) 



2m 



v' 2 + ^ M + y M cxt (r') 



+ E<W (n(r')-E ex t(r',f))^(r',f). (11) 

Correspondingly the equations of motion for the conju- 
gate momentum II and the transversal vector potential 
A read 



(d t + ft(i)cV)ri(r', t) = -— V x (V' x A(r',i)) 



(12) 



-n(r',t) + -P(r',t). (13) 



In eq. (|13|l we have introduced the transversal polar- 
ization P(r,t) = £ /li „#] l ( r ' < ) d /u'M r >*)- K is im- 
mediately obvious that the transformation to the ro- 
tating frame just amounts to the replacement dt — > 
dt + f2(t) 8^. For notational simplicity we will omit in the 
following the prime that indicates rotating coordinates. 

As we work in the Coulomb gauge we have II(r, f) = 
-D(r,i) |2^. Using this and D(r) = e E(r) + P(r) we 
find for the wave equation of the electric field in the ro- 
tating frame 



E(r,t) 



= -[d t 



n(t)d v P(r,t). (14) 



We 
the 



now introduce slowly varying variables for 
transverse field as well as polarization by 

E(r,i) = £ (x,r±,t)e- i ^" t - k ^ + h.a. and 

P(r,t) = V {+ \x,r ± ,t)e- i(uJ " t - k ^ + h.a., where 
x = Rip is the arclength on the circle. Restricting 
ourselves to propagation along the periphery of the in- 
terferometer we find within the slowly varying envelope 
approximation and by neglecting terms 0{VLR/c) 



d t + cd x 



S {+ \x,t) = - i ^V {+ \x,t). (15) 



The term proportional to the rotation rate fl is responsi- 
ble for the rotation induced Sagnac phase shift in the pure 
light case, i. e. without any influence from the medium 
polarization. As shown in [l4j and in the next section 
the polarization leads to an additional phase shift. 

Introducing also slowly varying amplitudes for the 
matter fields *i = <f>i, * 2 = $ 2 e~^ Upt ~ kpX ^ and # 3 = 
$ 3 e -i(Aa)t-Ate) with Alj = lj p - uj c and Afc = k p - kl, 

where fcj is the wavevector projection of the control field 
onto the x-axis, we find 



2?i-7i(x)J*i=«j;*2, (16) 
[V 2 + h{A 2 - kpSIR - F 2 (x)))$2 = hTLp®! 

+ m c $ 3 , (17) 

v 3 - v 3 (x) + n(A 3 - v k p nR)) $ 3 = m* c & 2 (is) 



with 

h 2 d 2 

V t ,=ihd t + —^+ih(nR + r ]t ,v ICC )d x . (19) 
2m 

Here we have used the definitions A 2 = u> p — lo 2 — w rec 
and A 3 = Alo — lo 3 — r) 2 tu ICC for the one- and two-photon- 
detuning including the recoil shift (o> rcc = hk 2 /2m). 
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Vrcc = hkp/m is the single-photon recoil velocity. We 
have also introduced the dimensionless parameter rj — 
Ak/kp which describes the momentum transfer from the 
light fields to the atoms in state |3) as well as the ab- 
breviation rj^ = 5^.2 + TjS^s. Finally the definitions 

Slp. c = — d PiC • E^J. for the probe and control-field Rabi 
frequencies were applied. The shortened wave equation 
l(T5j l and the matter- wave field equations lJTfj )) -lJT5 )l are 
the basis of the following study of the sensitivity enhance- 
ment of the light-matter-wave hybrid Sagnac interferom- 
eter. 



III. SAGNAC PHASE SHIFT AND INFLUENCE 
OF EXTERNAL TRAPPING POTENTIALS 




FIG. 2: Set-up of a Sagnac interferometer with ring-shaped 
trap configuration supporting a superfluid ultra-cold gas 
(BEC) . The symmetric interferometer set-up allows for a dis- 
tinction of rotational from linear acceleration. 



In this section we will calculate the stationary Sagnac 
phase shift for the hybrid interferometer in the perturba- 
tive limit of low probe-light intensities. In particular we 
will analyze the effects of a trap potential which confines 
the atoms to certain regions in the direction of the in- 
terferometer path. For simplicity we assume a constant 
rotation rate, i.e. fi = 0, and consider the stationary 
state. All atoms are assumed to be initially, i. e. before 
applying any probe-field, in the internal state |1). This 



amounts to set $ 2 0) (t = 0) = ^'(t = 0) = 0. In the 
perturbative limit of small probe field intensity, $i is 
not changed by the atom-light interaction, i. e. it obeys 
the equation 

(^.+mRd x + (ex - Vi{x))^\x) = 0, (20) 

where t\ is the energy in the stationary state. Assum- 
ing |fi c | ^> IA2I, fep|f2|iJ, |V2(x)|/fi one finds in first order 
from eq. (|T7ll 



*(0), 



*W (a) = _^*(o> (a0| 



(21) 



which amounts to an adiabatic elimination of the excited 
state. Using this and eq. (|18|l we find 



1 



iru 2 



h 

2m 



+i (VLR + r]v rec )d x + 



e 3 ~ V 3 (x) 



^ ] n p (x) (22) 



<^ 0) / . \ 

= V j^y-p [iVxeA In $1 + k p flRj £l p (x) 



fir 



(0) Vkf)2 

+ if SIR + V v rcc - i-d x In $( 0) ) d x 



ft p (x), 



constant energy term proportional to e\ — £3 has been 

dropped. One recognizes that the fields and 
and thus the medium polarization in first order of per- 
turbation follow straight forwardly from the solution of 
eq.J23. 

We will now consider two cases. In the first case no 
confining potential for atoms in state |1) is assumed, 
which is equivalent to translational invariance on a ring. 
In the second case, discussed later, a trapping potential 
in the longitudinal direction x is taken into account. We 
will see that both cases lead to quite different results. 



A. Periodic boundary conditions in state |1) 

Let us consider the case that atoms in state |1) do 
not experience any confining potential in the x-direction. 
Since x is the coordinate along the periphery of the inter- 
ferometer, this amounts to considering a ring-trap con- 
figuration with periodic boundary conditions $^ (x + 

2t:R) — (x). The principle set-up is shown in 
Fig. J2J. With V\(x) = 0, eq. has the eigensolu- 
tions 

n 2 h 2 



$^(a;)=$oexp 



e n = nhQ 



2mR 2 



where <l>o is a constant. Taken as a continuous function 
of n, the spectrum e(n) is a parabola with minimum at 



mSlR 2 



(23) 



This is illustrated in Fig. [3] Taking into account that n 
must be a positive or negative integer, the state with the 
lowest energy corresponds to n = as long as |n m j n | < 
1/2, i.e. as long as the Sagnac phase shift per round trip 
is smaller than tt. 



where we have in addition assumed two-photon reso- 
nance, i. e. A3 = 0. In deriving the second equation, 
which is useful for later discussions, we have made use of 
eq. (|2U|) and assumed equal trapping potentials for the 
internal states V% — V3. Furthermore an unimportant 



1. Bose- Einstein condensate 

We now assume that only the lowest motional energy 
state in the internal state |1) is initially excited, e.g. a 
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Bose-Einstein condensate in the ring trap. In this case 
there is a uniform phase over the whole ring and we can 

set 4^ (x) = $o. It is important to note that the parti- 
cles in the ground state do not attain a rotational phase 
shift in this case. This yields with eq. (|22|l 



(24) 



i(SlR + r)v lcc )$>o 



d x fl p (x) 



1 



(i) 



and $^°' ) into the 



Substituting the expressions for <I> 2 
stationary, shortened wave equation, ea. l|15|) . for the ex- 
pectation value of the probe-field expressed in terms of 



cd x + ikpQR ) O p (a 



WW 



(25) 




FIG. 3: Parabolic spectrum e„ with minimum at n m in = 
—mQR 2 /h when taken as a function of the continuous pa- 
rameter n. 



where g = di 2 y / ujp/2heoF, dn being the dipole matrix 
element of the |1) <-> |2) transition and F the transversal 
cross-section of the probe-beam, we find 



hd 2 

{ccos 2 9 + (r)v Tec + QR) sin 2 6\d x - isin 2 9—^ 

2m 

= -ik p QR( cos 2 9 + 77 sin 2 9 ) fi p (x). 



n p (x) 

(26) 



Here we have introduced the mixing angle 9 through 
tan 2 9 = g 2 g/\Q c \ 2 , where g — ($J^o) is the density of 
atoms in state |1). Eq. I|2t)l) has a very intuitive inter- 
pretation. It describes the propagation of the probe-field 
with the group velocity 



c cos 2 9 + n v Tec sin 2 1 



(27) 



in the rotating frame. The propagation of light in an 
EIT medium is associated with the formation of a dark- 
state polariton, a superposition of electromagnetic and 
matter- wave components [T^ . If we neglect the atomic 
motion, the group velocity of this quasi-particle is pro- 
portional to the square of the weight factor cos 9 of 
the electromagnetic part of the polariton. However, if 
the coherence transfer from light to atoms is accompa- 
nied by a finite momentum transfer of r]mv Tec , there is 
also a matter-wave contribution to the total group ve- 
locity (|27|l . This contribution is again proportional to 
the square of the weight factor sin 9 of the matter- wave 
part. Due to the admixture of the matter wave, the equa- 
tion of motion l|26|l attains a term corresponding to the 
kinetic energy of this component which leads to a dis- 
persive spreading of the probe field along its propaga- 
tion direction. This term becomes important in the limit 
tan 2 9 > tan 2 9 CI i t = c/v ICC , i.e. when the light wave es- 
sentially turns into a propagating spin polarization. The 
right hand side of eq. I|26l) describes the light and matter- 
wave contributions to the rotational phase shift. The 
matter-wave contribution to the phase shift is non-zero 
only if there is a finite momentum transfer, i.e. if r\ 7^ 0. 



In the limit of small rotation, |0|i? -C u gr , which is the 
case of interest here, eq. <|2(j|) can easily be solved. Ne- 
glecting the second-order derivative the equation reduces 
to eq. (11) of ref. 



where 



.27rflR / £(x 



mc r\ 



Ac + v ft^p £( x ) + r i 



cot z 6» 
cot 2 crit 



(28) 
(29) 



The last approximate equation is only valid for v gT <C c. 
When £ is large the group velocity is much larger than 
the recoil velocity, while £ approaching zero means that 
the group velocity is comparable to the recoil velocity. 
Eq. 128|) describes a phase shift of the probe field in a 
medium without absorption, which is canceled due to 
EIT. Hence two counter-propagating probe fields will ex- 
perience the Sagnac phase shift 



Acf> 



2kVLR 
Ac 



dx 



a*) 



QR 



£(jj)+?7 H/m J £(x)+i] 



dx ■ 



(30) 



This is the result obtained in [Tj]. It has two terms, a 
light-contribution and, if rj =/= 0, a matter- wave contribu- 
tion. Its most important consequence is that if the group 
velocity becomes comparable to the recoil velocity, i.e. 
for £ — > 0, the slow-light Sagnac phase approaches the 
matter- wave value! 



2. thermal gas 

The ground-state solution <lp(x) = 3>i°„ =0 = $0 = 
const, means that the atoms do not follow the motion 
of the trap. This is strictly speaking only possible if 
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the gas is superfluid. In a normal gas collisions with wall 
roughnesses and between atoms, which are not taken into 
account here, would accelerate the vapor particles in the 
initial phase of rotation. Eventually an equilibrium state 
would be reached where the atoms co-rotate with the 
trap. This can also be seen from a different argument. 
In a thermal state with ksT 3> Ml + h 2 /2mR 2 many 
states in the spectrum of Fig. [3| will be occupied. As a 
consequence the thermal gas in the ground state attains 
an average rotational phase (x — 2itR) 



A<p = -2n{n) 



2.7T 72>min — 



2-kVLR 2 
h/m 



(31) 



This is just the matter- wave Sagnac phase and is in sharp 
contrast to the case of a Bose-Einstein condensate, where 
the ground state does not acquire any rotational phase. 
Since now both, the ground state |1) and the excited 
state 1 2) attain the same Sagnac phase shift, the matter- 
wave contribution to the polarization is exactly cancelled. 
Thus the extention to thermal gases made in [14( is not 
correct. 

The need for a superfluid gas (e.g. BEC) in a ring trap 
puts restrictions to the achievable interferometer area. 
Although recently there has been substantial progress in 
realizing ring traps for BEC |18l 1 1 9| , the area achieved 
is only on the order of 1CP 1 cm 2 , which cannot compete 
with the values reached in fiber-optical gyroscopes. 



B. Effect of longitudinal confinement 

Let us now discuss the case of a longitudinal trapping 
potential for atoms in state |1), i.e. V\{x) ^ in eq. (J2DJ. 
In this case the substitution 



leads to the steady-state equation 

'h 2 dl 



n 2 R 2 + e 1 -V 1 (x) } f(x) = 



2m 



(32) 



(33) 



If one disregards the small centrifugal energy shift pro- 
portional to f2 2 , this equation is just the stationary 
Schrodinger equation for a particle in the trap potential 
V\ . The solution of this equation is independent of the 
rotation rate f2. (The inclusion of the centrifugal term 
would lead to a higher order contribution to the Sagnac 
phase, which we are not interested in.) If we substitute 
(I32JI into the second equation of Q22JI. one recognizes that 
all terms containing the rotation rate f2 vanish exactly: 



H 1} = ~ i ^£M (vvr ec (d x lnf(x))n p (x)-i^fl p (x) 



.&i 0) (s) 

1 \n r \ 2 



vv rcc - i— d x lnf(x) d x fl p (x). (34) 
m ' 



Substituting this into the shortened wave equation for fl p 
yields 



ccos 2 + (rjVrec - i—d x lnf(x))sm 2 0d x rtp(x) 
\ m / J 

hd 2 

isin 2 0^-%(x) = -ik p QRcos 2 0n p (x) (35) 



+ rjVrec sin 2 0^d x hi f(x)jQ, p (x). 



Neglecting the term with the second order derivative as 
well as those containing d x f(x), which amounts to as- 
sume that f(x) is a slowly- varying ground-state wave 
function of a sufficiently smooth potential, ea. l|35(l re- 
duces to 



d x In Q p (x) = — i 



2vrOi? 



cos 2 



Ac cos 2 + rj Vrcc sin 2 
2nnR 1 



Ac l+Tic/f 



(36) 



It is immediately obvious that only the light part of the 
Sagnac phase survives. Thus in the EIT hybrid gyro- 
scope a matter-wave contribution to the Sagnac phase 
only emerges in the absence of a confining potential or 
if periodic boundary conditions apply as e.g. in a ring 
trap. 

The physical interpretation of this result is straight for- 
ward. In the presence of a confining potential the atoms 
in state |1) are bound to the motion of the confining po- 
tential. Hence they acquire a rotational phase shift by 
following the motion of the trap attached to the rotating 
frame 17] . Atoms in state |2) acquire the same phase 
shift since they are in the same frame. Therefore, the 
polarization attains no Sagnac phase as it is a sesqilinear 
function of the wave- functions of states |1) and |2). This 
is in contrast to a superfluid BEC in a ring trap, where 
the order parameter does not acquire any phase due to 
the periodic boundary conditions as long as the rotation 
is sufficiently slow. 



IV. QUANTUM LIMITED SENSITIVITY OF 
THE SLOW-LIGHT GYROSCOPE 

We now want to calculate the sensitivity of the slow- 
light Sagnac interferometer in the case of periodic bound- 
ary conditions, i.e. in the absence of any confining po- 
tential in the propagation direction. For simplicity we 
consider the case r] — 1, i. e. perpendicular propagation 
directions of probe and control field. 

To determine the sensitivity we assume that the error 
in determining the Sagnac phase is entirely given by shot- 
noise quantum fluctuations. If coherent laser light or 
Poissonian particle sources are used the shot-noise limit 
of the phase measurement is given by 



Acj) n 



1 



(37) 
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where ?t,d = ^out^D is the total number of photons or 
atoms counted at the detector during the measurement 
time tv @. Here / ou t is the photon or atom flux. The 
assumption that the quantum noise limit is set by shot- 
noise is justified by two observations: First of all, it is well 
known that using nonclassical light or sub-Poissonian 
particle sources does in general not lead to an improve- 
ment of the signal-to-noise ratio in interferometry since 
at the optimum operation point the amplitude reduction 
due to losses is typically of order e _1 and thus quite sub- 
stantial. These losses tend to quickly destroy the fragile 
nonclassical and sub-Poissonian properties. Secondly, as 
has been shown in j2|| |2]j, atomic noise contributions 
in EIT-type interferometer set-ups are small and can be 
neglected. 

In the weak-signal limit discussed in the previous sec- 
tion, the Sagnac phase accumulated is independent of 
the signal field strength [T^ |. hence the signal-to-noise 
ratio could become arbitrarily large when the input-laser 
power is increased. In reality the Sagnac phase ap- 
proaches a maximum value at a certain optimum probe- 
laser power and decreases for larger intensities. The 



where 72 = 71+73- In the following we determine 
the Sagnac phase shift for arbitrary probe-field Rabi fre- 
quency based on the above set of equations and the short- 
ened wave equation. To derive a transparent expression 
for the rotationally induced phase shift further simplifi- 
cations are however necessary. 



optimum intensity is reached when the number density 
of photons in the EIT medium approaches that of the 
atoms. Thus in order to calculate the maximum sensi- 
tivity and to find optimum operation conditions we have 
to calculate the Sagnac phase to all orders of the sig- 
nal Rabi frequency £l p . Since in higher order pertur- 
bation the excited state |3) attains a finite population, 
decay out of the excited state needs to be taken into 
account. The decay leads to a population redistribu- 
tion among the states of the A system, see Fig. 0] It 
can also lead to loss out of the system. We will disre- 
gard the latter process however. Furthermore, we assume 
that the density of the considered medium is low enough 
that it is sufneent to describe the system by a set of 
equations for the single-particle density matrix elements 
t) = (®t{ x 'it)®i X (x,t)). Here p, v <E {1,2,3} 
denote the internal states. Since the medium polariza- 
tion is determined by the local density-matrix element 
Pi2{x, x, t), i.e. x' = x, we consider only local quantities. 
For the density matrix elements diagonal in the internal 
states we find the equations of motion 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



A. Non-local terms 

One recognizes from eq. <|41[) and (|42l) that the local off- 
diagonal matrix elements are coupled to non-local quanti- 
ties of the form (x)(d x ^ l/ (x))). These terms cause the 
build-up of coherences between different internal states 
and different positions, which are zero in lowest-order 



I 

d t pn(x,t) = 71/922OM) - iftp(x,t)p 2 i(x,t) +itt p (x,t)pi 2 {x,t) +flRd x pn(x 1 t), 

d t P22(x,t) = -72P22OM) +i%(x,t)p 2 i(x,t) - in p (x, t)p 12 (x,t) 

+ ifl*(x,t)p 23 (x,t) - ifl c (x,t)p32(x,t) + (flR + v rcc )d x p22(x,t), 

dtPzs{x,t) = 73/022(2:, t) - itt* c (x,t)p 23 (x 1 t) + iQ c (x, t)p 32 (x,t) 
+ (flR + v ICC )d x p 33 (x, t). 

Likewise we find for the local coherences 

dtpu{x, t) = -(i(A 2 + flRkp) + l2 /2)p l2 {x, t) + i n*(x, t) Pl3 (x, t) 
- iQ,*(x,t)(p 2 2(x,t) - p n (x,t)) + (SIR + v ICC )d x pi 2 (x,t) 

d t p\ 3 {x,t) = -(i(A 3 + flRk p ) + ~fi 3 )p 13 (x,t) ~ ittp(x,t)p 23 (x,t) 

+ in c (x,t)p 12 (x,t) + {flR + v ICC )d x p 13 {x,t) +v rcc ($l(d x $> 1 )} , 

d t p23(x,t) = i(A 2 - A 3 ) - j2/'2)p2 3 (x 1 t) - iQ p (x,t)pi 3 (x,t) 

-ifl c (x,t)(p 33 (x,t) - p 22 {x,t)) + (FIR + v ICC )d x p 23 (x,t). 

I 
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perturbation. We now want to argue that these terms 
can be neglected. To this end we consider eq. (I22[l again 
disregarding second-order derivates and set V3 = 63 = 0. 
Hence we have 



$ 2 = 



k D VLR 



in. 



$iSl p - i 



la 



^ x ($iO p ). (44) 



Substituting this into the steady-state version of the 
equation of motion for $1, eq. (|16(l . remembering that 
there is no confining potential for atoms in state |1) in 
the propagation direction, we find 



d x <&i(x) = -is 



[k p nR - i(JlR + v Tec )(d x In «„)] ~ 



nE(i + s) + w r , 



(45) 



IVI^I 2 

the probe field picks up a Sagnac phase shift, we have 



where s — \Q p \ /\Q C \ is a saturation parameter. Since 



d r In £1 



p 



—i—kpCiR. 

c 



(46) 



With the help of this we finally arrive at 
kMR 



d^i = -i- 



1 — a- 



((OR) 2 ). (47) 



As a consequence the term v Tec (& 2 d x &i) in eq. (|41|l is of 
the order of 



-ifcpfii? ( 1 



P12 



(48) 



and is thus negligible as compared to 72P12/2. Using 
similar arguments one finds that the term v ICC ($ld x $i) 
in eq. (|42ll is of the order of 



Pl3- 



(49) 



Since ideally the ground-state coherence is long-lived, one 
has 713 — * 0. Hence neglecting this term is not as straight 
forward as for ea. l|48|) . However, adiabatically eliminat- 
ing the fast decaying optical coherence P12 in eq. (|41|l and 
substituting the resulting expression into the equation of 
motion of P13, eq. I|42[) . yields a term proportional to 
|f2 c | 2 72yOi3/2 which is much larger than kpflRpi^. Thus 
also the term V gT {^\d x ^i) can be safely neglected. As a 
result of this approximation the density matrix equations 
d3HJ-<E3J are self-contained and local. 



B. Perturbation theory with respect to 
characteristic length 

In the following we assume one- and two-photon res- 
onance, i.e. A2 = A3 = 0, and solve the above system 
of equations for the coherence of the |1) <-> 1 2) -transition 
in steady state to all orders in £l p . The density matrix 
equations l|38|) - (|43|l can, neglecting terms proportional to 
QRd x , be written in compact form as 



p(x,t)= [M(x)+v TCC T>d x )p{x,t), 



(50) 




FIG. 4: A configuration in which the Rabi frequency Q p drives 
the 1 «-» 2-transition and fl c the 3 <-> 2-transition (solid lines). 
Radiative decay from the excited level to |1) or to |3) goes as 
71 or 73 respectively (dashed lines). The dephasing rate of 
the 1 — 3 coherence is denoted by 713 . 



where M and D are 9x9 matrices. Even under station- 
ary conditions we are still left with a set of first order 
linear differential equations with space dependent coeffi- 
cient. Thus in order to find an analytic solution further 
approximations are needed. To this end, we make use 
of the fact that the off-diagonal density matrix elements 
are only slowly varying in space. Let L and T be their 
characteristic length and time scales. Normalizing time 
and space to these units by £ = x/L and r = t/T, eq.j5U|) 
reads 



d T p=[M 



(51) 



where typical matrix elements of M = MT read as 
k p RflT, with k p R, \flT\ > 1 and those of D = DT are 
of order unity. Since v Tec T/L is typically small compared 
to unity we can apply a perturbation expansion in this 
parameter. 

In zeroth order we disregard the term containing D. 
Hence in steady state we have to solve Mpss = with 
the constraint X^Pmm^) = n ( x )-> which reflects the con- 
servation of probability. 

In first order we find 



zL 1) w = (. 



= X - Vr, 



M-'no^p^ix) 



(52) 



Here M is a reduced 8x8 matrix obtained from M by 
incorporating the constraint X)uPw( x ) = n ( x ) an< ^ 
is the corresponding zeroth-order density matrix. The 
explicit expressions of all matrices and vectors can be 
obtained from (|38|l - (|43|l in a straight forward manner. 
They are however lengthy and will not be given here. 



C. Steady state Maxwell-Bloch equation 

To obtain the rotationally induced phase shift we ex- 
pand eq. up to first order in the angular velocity fi 
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and use the time-independent Maxwell equation 1151) in 
the rotating frame 



(cd x - ik p QR)n p (z) = ~ig 2 Np s 2 



(53) 



To determine p s 2 \ we furthermore neglected terms 0(7 2 3 ) 
and 7i3^p with n 6 N since we assume a long-lived coher- 
ence between the two lower states |1) and |3). In addition 
to this we made use of the EIT condition tl 2 » 71371 
and assumed for simplicity 71 = 73 = 7. 

With these assumptions we arrive at the following ex- 
pressions for the real and imaginary part of the suscepti- 
bility which determine the dispersion and absorption of 
the medium 



x '(n p )= f 3- 1 —[i 



g 2 N 



2\2 



(Qi + Ifipl") 



n. 



k p c 9 N \n 2 c + \n p \y 



with 



I^pI 2 ) 3 ' 



(54) 



(55) 



(56) 



The imaginary part of the complex susceptibility 
x(^p) — x' + l x" can be further simplified. One can 
easily see that the absorption constant is bounded from 
above by 



k P x" < 



7i3 9 2 N 

C ft 2 



713 
^rcc 



r 1 . 



In this limiting case the following equation arises 

713 

c 



<9 z lnf2 p (z) 



• tan 2 9 + ik p x'(tt p ) 



(57) 



(58) 



The first term in eq. (|58|l describes absorption losses 
due to the nonvanishing decay of the ground-state coher- 
ence, the second term the rotationally induced or Sagnac 
phase. Since the saturation of the absorption for increas- 
ing probe-field intensities is not taken into account, the 
losses are slightly overestimated. 



D. Quantum limit of gyroscope sensitivity 

Solving the shortened Maxwell equation (|58|l for the 
probe field with the all-order susceptibility, ea. (|54*jl . we 
can now determine the minimum detectable rotation rate 
^min of the slow-light gyroscope. This is done by max- 
imizing the signal-to-noise ratio (SNR) of the interfer- 
ometer with respect to the system parameters and set it 
equal to unity. The relative rotational phase shift of two 
polaritons propagating in opposite directions is given by 



A^ sig = dxk p X '(Q,n p (x)) - X '(-n,n p {x)) . (59) 



Using this and eq. (|54|l we find 



A0 sig = 



Ac 



dx 



1 

(1 + S(3=)) 2 

h/m J ~~ £(ar) + {1+ i [x))3 



dx 



(60) 



where s(x) = \tt p (x)\ 2 /Q 2 is the saturation parameter in- 
troduced before, and £(x), defined in eq. (|29(l . determines 
the character of the polariton. One recognizes that the 
matter- wave part of the signal phase - the second line of 
eq. (|60|l - decreases for increasing input probe intensity. 
The light part - first term in eq. (|60|l - approaches a con- 
stant value in this limit. At the same time the shot-noise 
phase error 



A0 n 



1 



(61) 



is inversely proportional to |f2 p (0)| exp (— kLm), where 
Lm is the length of the medium and the probe-field's 
source is located at x = 0. As a consequence of the 
different dependence of the signal and noise terms on 
the probe field strength, the signal-to-noise ratio SNR 
= A</> s i g /A(/> no i SG has the qualitative behavior shown in 
Fig. |SJ For very large laser fields the SNR becomes ar- 




FIG. 5: Schematic dependence of SNR on input probe-field 
Rabi frequency. The dash-dotted line indicates the contribu- 
tion of the matter-wave term, the dashed line that of the light 
term. The solid line is the sum of both contributions. 

bitrarily large, as the light contribution to the Sagnac 
phase becomes intensity independent and the shot-noise 
level decreases steadily. For small probe intensities the 
SNR has a local maximum due to the saturation of the 
matter- wave phase shift. As the matter- wave contribu- 
tion to the Sagnac phase is orders of magnitude larger 
than the light contribution, extremely large input inten- 
sities would be required to exceed the sensitivity at the 
first local maximum. (Note that Fig. [S] is not drawn to 
scale.) We thus consider only this first maximum when 
determining the quantum-limited sensitivity of the slow- 
light gyroscope. 
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Although it is rather straight forward to calculate nu- 
merically, on the basis of the above given equations the 
minimum detectable rotational phase shift, we are inter- 
ested here in an analytic estimate. For this we make 
some simplifying assumptions: First of all we consider 
the propagation of polaritons through a homogeneous 
medium. We furthermore ignore the space dependence 
of the functions and s(x) in the expression (|6UI) for 
the signal phase, which amounts to replacing |f2 p (x)| by 
its input value |fi p (0)| = |fi p |. As will be seen this only 
slightly overestimates the saturation of the signal at the 
optimum operation point. We also ignore the satura- 
tion of the probe field absorption, which again merely 
slightly overestimates the probe field losses at the op- 
eration point. Finally we only consider the dominant 
matter-wave contribution to the signal phase. Thus we 
have 



A0 s ig = 



flRL 



(1 + * 



h/m £(1 + s) 3 + 1' 



(62) 



In order to estimate the signal-to-noise ratio SNR = 
A0 s ig/A0 no i SO we now express the the shot- noise expres- 
sion (|61J) in terms of the parameters £ and s. The number 
of probe photons at the detector can be written in terms 
of the probe-field Rabi frequency at the source via 



p u t _ 2e Q Fc fnn p (o) 

hu>„ hui n 



te 



(63) 



where F is the cross-section of the signal beam, t the de- 
tection time interval, and k — 713 / (v Tec t;) the absorption 
coefficient introduced before. The radiative decay rate 
7 = 71 and the dipole matrix element |d p | contained in 
the Rabi frequency fi p (0) are related through 



(64) 




i. e. according to the Einstein A-coefficient [2f 
straight forward calculation we find 

n D = F gv Icc ttse- 2a/ Z 



After a 



(65) 



where g is the density of atoms in the EIT medium, and 



713 £ 



(66) 



characterizes the absorption due to a finite lifetime of the 
ground-state coherence. Since typical values of 713 are in 
the kHz regime and v rcc ~ 1 cm/s, a is typically large 
compared to unity for Lm ^> 10 -3 cm. With the above 
expressions we find for the signal-to-noise ratio 

SNR = —{Fgvr ec t y 



H/m 

£W/3(i +s ) /£ 

£(l + s) 3 + l 



(67) 



The first two factors in eq. (|67(l are the signal-to-noise 
ratio of a pure matter-wave gyroscope with interferom- 
eter area A = RLm and a flux of atoms corresponding 
to a density of atoms g passing through an area F with 
recoil velocity v rcc . In conventional atomic interferome- 
ters based on cold or ultra-cold atoms the flux that con- 
tributes to the interference signal of the device is rather 
low. It is on the order of 10 8 atoms/s in comparison with 
10 16 photons/s in a conventional fiber-optics gyroscope 
[Tlj . However, in the case studied here, the flux can be 
at least two orders of magnitude higher than in an atom 
interferometer. 

The second factor can be modified by optimizing the 
probe field strength (s) and the group velocity in the 
medium (£) . In Fig. we have plotted the optimum 
values of s and £ derived by maximizing the signal-to- 
noise ratio for different values of the loss parameter a. 



0.8 
0.6 

0.4 

0.2 



1 2 3 



FIG. 6: Optimum values of 



\n p (0)\ 2 /n 2 c and £ 



v gI /v lac — 1 for different values of the loss parameter a = 
Ji3L M /v TCC (= 0.05,0.5,5,50,500,5000). For large values of 
a the optimum values are s opt = 1/3 and £ op t = 2a. For small 
values of a there is only a small deviation in the optimum pa- 
rameters. 



One finds that in the typical parameter regime a 1 
the maximum SNR is attained for 



'opt 



and £ 



opt 



2a. 



(68) 



Note that this approximation is still quite good even 
when a = 1. The optimum group velocity is given by 



„°pt — 



2713 L 



2713-^M, 



(69) 



i.e. a maximum SNR is achieved if the velocity is chosen 
such that during the propagation over the entire medium 
length Lm, a fraction of 1/s/e of the initial polariton 
gets absorbed. Setting SNR = 1 we eventually obtain 
the minimum detectable rotation rate 



h/m 



1 



.4 



(F Q u rec t) 



1/2 



/V5 



(70) 



where / « 7.2 is a numerical prefactor resulting from the 
optimization of the term in the second line of ea. (|67|l . 
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Apart from the term ^[a and the unimportant numer- 
ical prefactor /, the minimum detectable rotation rate 
corresponds to that of a matter-wave interferometer 
with atoms propagating at recoil velocity. The densi- 
ties achievable in the present set-ups, e.g. if we con- 
sider BECs in ring trap configurations, are however much 
larger than those in typical atomic beams. 

To be more precise we give an estimate for the mini- 
mum detectable rotation rate of the slow-light gyroscope 
achievable with current technology. To this end we con- 
sider two state-of-the-art circular waveguides for Bosc- 
Einstein condensates [Til IT^|. Furthermore, we assume 
that the atomic density of the BECs is p = 10 14 cm~ 3 
with a cross-section (smaller circle of the toroidal BEC) 
of F « 1CU 2 cm 2 . In case of the work of S. Gupta et 
al. ^| the diameter of the larger circle of the toroidal 
waveguide is < fc» p ta ~ 3 mm and in the case of A. S. 
Arnold et al. it is ^Arnold ~ 96 mm. Hence, we find 
in the first case the minimum detectable rotation rate 
to be ^" n pta » 1.4 x 1(T 9 s- 1 Hz" 1 / 2 and in the latter 
case f^™°id rj 1.4 x 1(T 12 s" 1 Hz" 1 / 2 . These values 
should be compared to the state of the art which for op- 
tical gyroscopes is 2 x 10~ 10 rad s _1 Hz -1 / 2 [2£j and for 
matter-wave gyroscopes 6 x 1CU 10 rad s _1 Hz -1 / 2 [30|. 

V. CONCLUSION 

We have analyzed in detail a novel type of light-matter- 
wave hybrid Sagnac interferometer based on ultra-slow 
light in media with electromagnetic induced transparency 
(EIT) proposed by us in 0] . In particular the influence 



of confining potentials was investigated and the shot- 
noise limited sensitivity and the minimum detectable ro- 
tation rate determined. By combining features of light 
and matter- wave devices the hybrid interferometer yields 
a minimum detectable rotation rate which is potentially 
better than the currect state of the art by up to two or- 
ders of magnitude. We have shown that as opposed to 
claims in earlier proposals for slow-light gyroscopes [3l| . 
it is not sufficient to utilize only the dispersive properties 
of EIT-media to achieve an enhancement of the rotation 
sensitivity. It is rather necessary to employ simultane- 
ously coherence and momentum transfer in the associ- 
ated Raman transition. Moreover we have shown that 
the medium has to be prepared in a state in which it 
does not acquire any rotational phase shift. This can 
be achieved for example by using a superfluid BEC in a 
ring trap as EIT medium. The requirement for periodic 
boundary conditions reduces the potential of the hybrid 
interferometer idea as compared to the statements in |l4| 
as it is not possible to build large area interferometers un- 
der this condition with current technology. However, the 
potential large flux of the proposed slow-light interfer- 
ometer leads to a substantial reduction of the shot noise 
as compared to state-of-the-art matter-wave gyroscopes 
and thus leads nevertheless to a substantial sensitivity 
enhancement. 
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